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I. INTRODUCTION 



Entanglement has been extensively investigated in the recent past, as it is a critical re- 
source for quantum information processing [l|. One model of quantum computation, the 



one-way quantum computer 



2], relies explicitly on entanglement. The resource of entangle 




ment is not at all random pure quantum state is typically highly entangled 

fact there is so much entanglement in typical random pure states that recent studies 
find them not to be useful for one-way quantum computation. This motivates the question 
of studying subsets of states with a control over the amount of entanglement available. 

It is well known that most of the entanglement in many body quantum systems is mul- 
tipartite. In random pure states of iV qubits, we need to consider blocks whose total size is 
at least larger than N/2, for them to be entangled j?|. This being the case, entanglement in 
smaller blocks is nearly impossible to observe. Previous studies have shown how rare it is to 
have two qubits entangled in a many qubit random pure state 3, [8]. In this paper it is shown 
that there is a surprising connection between the number of up-spins or particles present in 
definite particle states and entanglement. Thus producing definite particle random states, as 
defined below, may allow control over the type of entanglement that is desired. For instance 
if two qubit entanglement is to be obtained, it is shown that typical three-particle states 
will render this nearly impossible to achieve. The border between probable and improbable 
is described by a transition from an algebraic to an exponential decay, which is typically 
obtained at phase transitions. Further, the approach presented in this paper might shed 
light on methods that are applicable to a wider class of problems in the area of quantum 
complex systems. 

Random pure states, or "full" random pure states, belong to the ensemble of states that 
are uniformly sampled from the Hilbert space, with the only constraint being normalization, 
in other words sampled from the unique Haar measure. Such states arise for instance in 

n n 

mesoscopic systems [9(, nuclear physics [10[ etc. and have been modeled as eigenfunctions 



of random matrices from the usual Gaussian ensembles Uj. There have been studies that 
explore how to efficiently produce operators with statistical properties of random matrices 
12| . Classically chaotic systems have long been known to exhibit such states in their quan- 
tum limit, and studies of entanglement in quantum chaotic systems often take recourse to 



random states 



13 
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The ensemble of interest in this work is taken to be the one where all the vectors that are 
constrained to be in the definite particle subspace are equally likely, subject again only to 
the constraint of normalization. This is equivalent to assuming that the Hamiltonian in the 
definite particle symmetry reduced subspaces are full random matrices. In this paper the 
random matrices are taken to be of the GOE (Gaussian Orthogonal Ensemble) H[ type and 
hence the states studied are real. Any study that uses random states or random matrices is 
justified as a baseline with which to compare realistic Hamiltonian systems that may involve 
interactions in a complex, nonintegrable manner. In particular it is possible that due to the 
few-body nature of typical interactions, ensembles such as the Embedded GOE [15[ will 
be of interest. However we believe that studying a usual ensemble like the GOE will form 
a baseline for entanglement studies of a rather large class of physically important systems 
which conserve particle number, or total spin. 

A definite particle state is a random pure state in a fixed S z subspace formed by the basis 
vectors of the Hilbert space, which, when expressed in the spin-z basis, have a fixed number, 
say /, of "ones" , or spin ups. Clearly many Hamiltonian systems including spin models such 



as the quantum spin-glass 



16| . or the disordered Heisenberg chains [l7] are potential places 



where such states can occur as eigenstates. The number of particles allows to add complexity 
to the states in a systematic manner, and interesting properties for entanglement unfold in 
the process. Studies of entanglement in two-electron systems for instance are found in 18]. 
One other class of problems where there is potential to see the kind of transitions noted in 



this paper is in the study of site-entanglement of fermions in a lattice [19j , where the total 
number of qubits of this paper will be translated to the total number of sites, the number of 
particles will be the number of fermions and the block will refer to the sites within which the 
entanglement is found. It may also be noted that translationally invariant definite particle 
states with highly entangled nearest neighbor were constructed as "entangled rings" 20]. 

A previous study of entanglement in random one-particle states showed that the averaged 
concurrence between any two qubits scales as 1/N 2l\ . Thus with increasing number of 
qubits entanglement between any two still remains considerable, although decreasing, in 
contrast to a full random state. In this paper it is shown that for random two-particle states 
the average entanglement between two qubits scales as 1/N 2 , while for three-particle states 
this becomes exponentially small, as it goes as exp(— Nln(N)). Thus when the number of 
particles exceeds two a transition is seen in the entanglement between two qubits. It maybe 
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noted that for full random states it is not precisely known how such an entanglement scales 
with the number of qubits. 

It is possible to generalize the results of concurrence between two qubits to entanglement 
within the block A having m qubits of the system for instance by studying the log-negativity 
measure |22| . Numerical and some analytical evidence points to the plausible result that the 
entanglement decays with N algebraically if the number of particles in the subspace (I) is 
less than or equal to the block- length (m). Once again the decay of entanglement becomes 
exponential when the number of particles exceeds the block-length. A study of the density 
of states of the reduced density matrix also shows a transition when the number of particles 
exceeds the block size; namely a divergence at the zero eigenvalue is replaced by a vanishing 
density. This paper studies this divergence and how this impacts the partial transpose in 
such a way that entanglement undergoes the kind of transition that is discussed herein. 

The structure of the paper is as follows. In section [Til details of definite particle states 
and reduced density matrices of blocks of qubits are given. In section II III the entanglement 
between two qubits is studied as a function of the number of particles present in the state 
and as a function of the total number of qubits. Detailed analytical results are derived 
utilizing concurrence as a measure of entanglement. To demonstrate that the transition is 
independent of the particular measure of entanglement, as well as to facilitate generalization 
to larger blocks, the log- negativity between two qubits is also considered here. In section 
JV] entanglement among larger sets of qubits is studied by means of log- negativity and a few 
analytical and several numerical results are presented. In section [V] the density of states of 
the density matrix is studied and it is shown that there is a transition in its character as 
the number of particles exceeds the block size. Evidence is presented that this results in, or 
is reflected as, the entanglement transition. 



II. DEFINITE PARTICLE STATES 

A definite /-particle state is best written by grouping states with a given number of 
particles present in one block, say A, and its complementary block, say B. Let the number 
of qubits in block A be m and let I > m. Label the states by the number of particles (or 
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total spin S z ) in subsets A and B to write: 



f N — m\ ( m\ 

\ 1-1 



(T) 

E 45 ) ioui%+ E E^i 1 )^- 1 ^ 

i=i i=i »=i 



(N-m\ 
\ l-m J 



+ ...+ E c™\m) A \l-my B . (1) 
i=l 

The reduced density matrix of the subsystem A denoted pa, is the state of the block of 
qubits we are interested in studying. These blocks correspond to having a given number 
of particles, k, in the subsystem A and can be identified with one of the terms in the 
expression for the state. Further, each of these blocks can be written as Gk = QkQ\ where 

(k) 

Qk is a matrix whose entries are the coefficients c\j of the state. The dimensions ( Throws x 
^columns) of, the in general rectangular matrices, Qk and the square matrices Gk are given 
by 

*■*-(:) KIT)'— (:)*(:)■ < 2 > 

The condition that the trace of a density matrix is unity implies that ^ k Ti{QkQ\) = 1. To 
construct the ensemble of /-particle states, draw all the M = (^) coefficients a-j from the 
normal distribution N(0, 1) and normalize them so that the trace condition is met. This is 
equivalent to choosing them uniformly with the only constraint being normalization 

The case when the number of particles I is less than the block length m can be similarly 
written: 



23|. 



^T) (T) 



w= E c ?j\ q )a\ 1 )b+ E E^i 1 )^- 1 )* 

j = l j = l 8=1 

(T) 

+ --- + E c ?l^|0) B . (3) 

In this case the last non-zero block has dimension ('?), however it has only one non-zero 
eigenvalue, as the corresponding density matrix QiQ\ is that of a unnormalized pure state. 
The sum of the dimensionality of the Gk block matrices in this case is less than 2 m and 
hence there are exact zero eigenvalues, in other words, the density matrix is rank-deficient. 
For example in the simple case of one-particle states, the two qubit density matrix has one 
exact zero eigenvalue. In fact it is easy to see that in general for one-particle states, the 
number of non-zero eigenvalues of the reduced density matrix of any number of qubits is at 



6 



most two. From the dimension of the Qk matrices above we can formally write the minimum 
number of exact zero eigenvalues of the reduced density matrix for the case when I < m. 
Written in terms of rank: 



In the case of I > m the rank of typical reduced density matrices is full, that is the rank 
is 2 m . This is seen by examining the dimensions of each of the blocks Gk- From Eq. (j2J) it 
follows that typical Gk will not have exact zero eigenvalues if 



A straightforward but nontrivial calculation shows that this is always the case if N > 21. 
The equality is possible only when I = m. It is indeed assumed that N > 21 throughout this 
paper without any loss of generality, due to the particle-hole symmetry. In the special case 
of N = 21 and I = m all of the Qk matrices are of square type and the Gk are formed then 
from symmetric states. This has implications for the density of states as will be discussed 
later in this paper (see section [V}. With these details about the structure of the reduced 
density matrices of definite particle states, we now turn to a study of entanglement in them. 

III. ENTANGLEMENT BETWEEN TWO QUBITS IN A DEFINITE PARTICLE 



The reduced density matrix of a block A with m = 2 qubits in a pure state of I > 2 
particles with the total number of qubits equal to N can be written as: 






(5) 



STATE 



/ 



a 00 



an ai2 







(6) 



Pa 



a* 2 a 2 2 
\ 



033/ 







where, a 00 = Ef^C) 2 , a 33 = ££i(c<?) 2 , and 
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Here fa = ( ,_ 2 ), i = 0,1, 2. The results presented in this work deal with real coefficients, 
c ij = Hi i a situation that would be relevant for example for systems with time reversal 
symmetry. The central features, including the scaling, remain the same in the complex 
case. Also note that while the above expressions have been written when I > 2, it is 
straightforward to write the same when / = 1, the case of 1-particle states. As a previous 



work has dealt with 1-particle states 



2l|, this is not considered further, however a detailed 



analysis of log-negativity in this case is presented later on in this paper. 



A. Concurrence between two qubits 



Concurrence |24j is a measure of entanglement present between two qubits such as those 
in the subsystem A. The above structure in Eq. ([6]), greatly simplifies the expression for 
concurrence [25 ] 

C = 2 max(|ai 2 | - V a ooa33, 0), (8) 
and this allows for analytical estimates to be made, in contrast to the case of a full random 

(k) 

state. Due to the large number of coefficients q- involved, it is a good approximation 
to assume that the normalization constraint is only important to set their scale and that 
they are otherwise independent. This implies that these are i.i.d. random variables drawn 
from the normal distribution N(0, 1/JV); recall that J\f = (^). The random numbers are 
ingredients for the random variables aoo, cti2 and 033 that determine the concurrence and 
hence the entanglement. 

The approach to finding the mean concurrence will be to first estimate the probability 
that it will be nonzero. The term a\2 involves a correlation between two strings of normally 
distributed numbers, each of length fa, the two rows of the matrix Q\ in Eq. ((7j). On the 
other hand aoo maybe taken to be effectively its average and considered to be non-fluctuating, 
as it is a sum over ~ N l random terms. The following approximation then ensues: 

Pr(C>0)^Pr(|a 12 |- v ^oT v ^>0). (9) 

The distribution of | ai2 [ , P12, is of central importance and can obtained from, for example, 
from the probability density function of one of the marginals of the Wishart distribution for 
correlation matrices 



261 ] . Suppressing the calculation, the result is 

K v {Mx) ( N_ 
2 



P 12 (\a 12 \=x) = 2Af J}. > — , (10) 
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where K v (x) is the modified Bessel function of the second kind, and v = (/ii — l)/2. The 
distribution of 1/033, -P33, follows from that of (square root of) a chi-squared distribution 
with fi 2 degrees of freedom: 

A/W2 

P 3 3(^= y )= 2w / lr(/t2/2) y »-'e-^/-. (11) 
Thus in view of the approximation above it follows that: 

POO POO 

Pr(C > 0) = / P 33 {y) / P 12 (x) dxdy. (12) 

J0 J y/ {a Q0 }y 



This can be evaluated by steps which are outlined here: (a) change variable y to ^2y/Af, 
and Nx to 2; (b) use the integral representation K u (x) = J °° e~ xcosht cosh(z4) dt, and change 
variable from xcosht to x. This leads to the exact expression (given the approximation in 
Eq. ®) 

Pr(C > 0) = (3 [°° COS ^l\ f" y*- x e^T{ V + \,jTvca&t)d V dt. (13) 
Jo cosh t Jo 

Here = 2~ u+l /^T(fi 1 /2)T(fi 2 /2), and 7 = 2(a 00 )Af. While further simplification is 
possible, for example by expanding e~ y , it is expedient to seek a non-trivial upper bound 
that reveals the nature of the decay with N, the number of qubits. A careful examination 
of the integrands indicate that this can be most easily achieved by using e~ y < 1 and thus 
effectively removing the exponential from the integral. The remaining integrals can be done 
exactly to give the first inequality below, while the second follows from standard inequalities 
for ratios of gamma functions: 

Pr < c > °> < r ,/ 2 0F r (f ) r(f + 1) * ^ \ 7 ) /5— T (14) 

where ?y = 1 + (/i2 — 2)/ (2/^). Note that when the number of particles is much less than the 
number of qubits, r] w 1. For the case of two particle states, 1 = 2, the first inequality yields 

Pr(C > 0) < —^=, (15) 



as /ii = iV — 2 and fi 2 = 1- The inequality is valid for large AT, especially as the value of 
(oqo) is taken to be 1. As a matter of fact that this is an excellent estimate by itself is seen 
in Fig. (pp. 
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log 2 N 

FIG. 1. The scaling of Pr(C > 0) for random two particle states with iV qubits. The dashed line 
is of slope —1/2, the circles are from numerical simulations, while the solid line is the estimate in 
Eq. (fT5jh Inset shows the average concurrence, the dashed line is of slope —2, the dashed-dot line 
is the upper-bound while the solid line is the estimate in Eq. (117p . 



The two particle case is of special interest and can be essentially derived from simpler 
formulae, if it is observed that the fluctuations in 033, arising from a single realization of the 
random variables, is more than the others. Note that: aoo ~ ^o/M ~ 1, ^33 ~ /J>z/J^f ~ t/N 2 , 
and 1 012] 2 ~ fii/Af 2 ~ 4/iV 3 . Hence typically the concurrence will indeed be zero. Replacing 
the average values for the fluctuating aoo an d 1 0,12 1 results in Pr(C > 0) ~ Pi(y/a^ < 
I72) — coinciding with the upper bound just derived. The average value of |a.i2 1 

is used, rather than the (square root of the) average of |cti2| 2 ; the exact distribution can 
be used to show that (|ai 2 | 2 ) = f ( ( 1 012 1 ) 2 ) - Thus for two particle states the probability of 
concurrence being positive decreases algebraically, in contrast to the one-particle case when 
P(C > 0) = 1, as a 33 = 0. 

For I = 3, three particle states, a completely different behavior is obtained as fix ~ N 2 /2, 




10 



fi2 ~ N, and 7 ~ N 3 /3 which result in 

Pr(C > 0) < ^/^exp (-^ log(JV/3)) . (16) 

Unlike the two-particle case the probability that the concurrence is positive decreases at 
least exponentially with the number of qubits, see Fig. fl2]). Another new feature is that it is 
quite essential to take into account the fluctuations in both | ai2 1 and in 033. Ignoring say the 
fluctuations in 033 results in much smaller estimates of the probability than what is found. 

When I > 2, but still much less than N, the upper-bound in Eq. (THI) does not estimate 
the probability accurately. While it can be made tighter, this is indeed a good bound as it 
is simple, decreases with N, and shows the advertized transition in the entanglement as one 
particle is added to a two particle state. It will be seen that the entanglement hitherto shared 
between two qubits will now be available for three-body and multi-party entanglement. 

If the fraction of particles p = l/N is of order 1 (and less than 1/2), the states are 
"macroscopically" occupied; employing the approximation that ~ e SN where S = 

— pln(p) — (1 — p) ln(l — p) is the binary entropy corresponding to probability p, results in 
the upper bound Pr(C > 0) < die~ SN ^ 2 e~ d2eSN , where d\ and di are positive constants of 
order 1. However the upper-bound in Eq. ffl4|) has to be used with caution as it can be 
rendered trivial if (2fiii]/ / -f) > 1, and consequently 62 becomes negative. Thus for N = 10 
qubits and I = 5 particles the upper-bound ~ 2.2 is trivial while for / = 4 it is ~ 1.5 x 10~ 5 . 
While N — 12, I — 6 results in a trivial bound, I = 5 results in Pr(C > 0) < 3.3 x 10~ 15 . 
Similarly when N = 14 and / = 6, the upper-bound is ~ 1.6 x 10~ 43 , it is improbable that 
two qubits will be entangled. 

The mean concurrence, E(C) is now estimated. In the two particle case for instance 

E(C) ~ 2<|o la |>Pr(C > 0) ~ - JjJL (17) 

A more general estimate is possible as E(C) = E[2(x — a/ (a 00 )y)Q(x — \J (a 00 )y)]. Using the 
distribution Pi2{x)P^{y) and following the same steps as outlined for the probability above 
it follows that 

where rf — 1 + (/i2 — l)/(2/^i). In the two particle case this gives E(C) < 8/y/iiN 2 , which 
is quite close to the estimate above. The exponential decay for three or more particles is 
manifest. The mean concurrences are shown in the insets of Figs. ([I]),©. 
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6 8 10 12 14 16 
N 



4 6 8 10 12 14 16 

N 

FIG. 2. The probability Pr(C > 0) for three particle states, and the average concurrence as 
number of qubits N is changed. Note that the y-axes are on a logarithmic scale. The circles are 
from numerical simulations while the solid line in the case of Pr(C > 0) is from an exact numerical 
evaluation of Eq. (fT2|) . 



B. Log-negativity among two qubits 



The vanishingly small two qubit entanglement for more than two-particle states (I = 2) 
goes into multiparty entanglement. A measure of entanglement that can be easily extended 
to a subsystem having more than two qubits is the log-negativity 22J and is given by 
Et,n{p ab ) = log (| |p^ B | |), where ||p r || is the trace norm of the partial transpose matrix 



27| . When log-negativity is zero the state is said to have positive partial transpose 



(PPT) and in that case it is either separable or bound entangled 28[|. When log- negativity 
is greater than zero the state is said to have negative partial transpose (NPT) and in that 
case it is entangled. 

On studying entanglement in a block length of 2 we get the entanglement between two 
qubits as measured by log-negativity. This decays algebraically as 1/iV 3 ' 5 in contrast to the 
1/N 2 behavior of the concurrence for the case of two particles, but becomes exponential 
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when the particle number is increased to three or more. See Fig. ([3]) for details. Thus on 
using a different measure of entanglement while indeed the exponents change the qualitative 
nature of the decay with the number of particles remains intact. It is also useful to contrast 
the case of 1-particle states and therefore log-negativity is now derived between any two 
qubits for both 1- and 2-particle states. 



1. Block of 2 qubits and 1-particle states 



In this case the reduced density matrix is block diagonal consisting of two square block 
and is given as follows: 



Pa 



where, a 00 = E*=i 2 ( c i?) 2 > and 



an a\2 



fa 00 0\ 

an &12 

a* 12 a 22 

\ oj 



QiQ\, and Qi 



"12 



a 22 



„Ci) 
ii 

-21 



where to remind the reader the coefficients c A A £1X6 £IS defined in Eq. (JTJ 
Partial transpose (PT) on the second qubit of pa results in 

f a 00 a 12 \ 



Pa 



(19) 



(20) 



(21) 



an 
a 22 
\a* 12 / 

The eigenvalues of p v A (as always in this paper, for the case of real coefficients c|^) are 
A± = ± v/(a o) 2 + 4a? 2 )/2, an and a 22 . The only negative eigenvalue is A_. From 
the assumptions of randomness of the state, for one-particle states we see that aoo ~ 1, 



'12 



1/N 2 . Using this the negative eigenvalue can be approximated by — ai 2 /aoo 



t 12- 



The log-negativity is given by Eln = log(l — 2 J2i u i) ~ — 2 J2i u i ~ ^ a i2 ~ 2/iV 2 , where the 
sum is over all the negative eigenvalues (cuj) of p T A and the approximation holds good since 
the tUj's are much smaller than 1. Indeed one finds that this estimate is in good agreement 
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with the numerical results shown in Fig. fl3]). Note that the average concurrence between 
any two qubits for 1-particle states scales as 1/N |2l| . 



2. Block of 2 qubits and 2-particle states 

In this case reduced density matrix is block diagonal consisting of three square block and 
is given as in Eq. fl6]). Partial transpose on the second qubit of pA results in 

(a m a 12 \ 

an 

a 22 

\a* 12 a 33 / 



P T a 



(22) 



The eigenvalues of this are 
A± = 



aoo + «33 ± V (aoo + «33) 2 - 4(a oa 33 - af 2 ) 



and the pair of a n and a 22 . Again the only possible negative eigenvalue is A_ which occurs 
when aoo&33 — of 2 is negative. In the case of two-particle sates, as we see earlier, aoo ~ 1 5 
4/A 3 and a 33 ~ 1/N 2 and thus the eigenvalue A_ can be approximated by (aoo&33 — 

a\ 2 . Thus E LN = 2a\ 2 Pr((a oa33 - al 2 ) < 0). Using Eq. (1I5|) 

16y/2 1 



a 12 ~ 

a 12)/( a 00 + CI33) 
we find that 



^00^33 



E 



LN 



7T iV 3 - 5 

which is in good agreement with the numerical results shown in Fig. ([3]). In the case of a block 
of 2 qubits and 3-particle states one finds that the log-negativity scales exponentially with 
the total number of qubits iV as shown in Fig. ([3]). This follows on using the exponentially 
small probability for the concurrence to be positive (see Eq. ( fl6l) ) and a similar analysis as 
above. 



IV. ENTANGLEMENT AMONG LARGER SUBSETS OF QUBITS 



While the previous section has dealt exclusively with a "block" of two qubits, here we 
take larger subsets of qubits to belong to block A. The Log-negativity measure will be used 
once again. A transition similar to the one above is exhibited for the entanglement between 
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FIG. 3. Scaling of log-negativity {Eln) m a block of 2 qubits, with the total number N of qubits 
for (left) one, two-particle and (right) three-particle states. In figure on left dashed line, dashed-dot 
line are having slope of -2, -3.5 respectively. 

a qubit and the other pair when a block of 3 qubits is considered. Algebraic decay of the 
log-negativity for I < 3 is replaced by exponential decay for I > 3, see Fig. The decay 
with particle number is algebraic and the exponent is the "slope" in Table HI Further results 
for block lengths of 4 are presented in Table HIl in which case there are two distinct types of 
partitions, entanglement between two pairs of qubits (denoted as 2+2) and between a triple 
and a lone qubit (denoted 3+1). However the numerics becomes considerably more difficult 
thereon, and the slopes given may not be entirely converged. However, the transition from 
algebraic to exponential is a robust feature. 



TABLE I. Block length 3. 



Particle # (I) 


Decay with number of qubits N. 


1 


Alg.: slope = -2 


2 


Alg.: slope = -3 


3 


Alg.: slope = -4.5 


4 


Exponential 



It is possible to extend the analysis of the log-negativity of two qubits to the case of a 
block of three qubits and derive some of the exponents as stated in the table. Now large 
N formulae for log-negativity in the case of a block of 3 qubits in 1-particle and 2-particle 
states are derived. It is shown that the log- negativity decays as 4/iV 2 , 16/iV 3 for these cases 
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TABLE II. Block length 4. Cases of (2+2,3+1) 



Particle # (I) 


Decay with number of qubits N. 


1 


Al? • slones —('-2 1 -2 11 


2 


Alg.: slopes = (-2.1, -3.1) 


3 


Alg.: slopes = (-4.1, -4.1) 


4 


Alg.: slopes = (-5.7, -5.7) 


5 


Exponential 




FIG. 4. Scaling of log-negativity (Eln) m a block of 3 qubits, with the total number N of qubits 
for (left) one, two and three-particle and (right) four-particle states. In figure on left dashed line, 
dashed-dot line and simple straight line are having slope of -2, -3, -4.5 respectively. 

respectively. 



1. Block of 3 qubits and 1 -particle case. 

In this case the reduced density matrix is block diagonal consisting of two square blocks, 
one of them (a o) being just a number: 



Pa 





0lx3 


0lx4 


3 xl 


QiQ\ 


3 X3 


04x1 


04x3 


04x4 



(23) 
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where a 00 = ^f =1 3 (cii) 2 , O pxq are zero matrices with dimensions p x q and 



QiQx 



( a\\ CL12 Ol3 ' 

a l2 a 22 a 23 
\ a 13 a *23 fl 33 J 



Ql 



c ll 

c (1) 

c 21 

v c (1) / 

V C31 / 



(24) 



Being 1-particle states these have only two nonzero eigenvalues in general. On PT it is seen 
that there are four nonzero eigenvalues. Partial transpose on the third qubit of pa results 
in 



P T A 



^00 











ai2 


ai3 





o\ 





a n 


























Ct22 


Ct23 




















a 23 


«33 














a *12 























a *13 



































































( V 



(25) 



The nonzero eigenvalues of p r are an, 022 + 033 and 



A, 



a 00 ±Ja 2 00 + 4{a 2 12 + a 2 13 ) 



Note that there are correlations in the entries of the QiQ\ matrices, such as 022033—^23 = as 
the state from the density matrix is constructed is an (unnormalized) pure state. Here aoo ~ 
1, a\ 2 an d a i3 ~ 1/N 2 . It can be seen that only one of the four nonzero eigenvalues is negative 
and it is A_. Using this, the negative eigenvalue can be approximated as — (af 2 + ai3)/aoo ~ 



(0^2+013)- The log-negativity is therefore given by E 



LN 



This estimate is in good agreement with numerical results as shown in Fig. 
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2. Block of 3 qubits and 2-particle case. 



In this case the reduced density matrix is block diagonal consisting of three square blocks 
and is given as follows: 



Pa = 



^ o-oo 


0lx3 


0lx3 


Olxl\ 


03x1 


QiQ\ 


3 x3 


3 xl 


3 xl 


3 x3 




3 xl 


\0ixl 


0lx3 


0lx3 


OixJ 



(26) 



where a 00 = EZii^) 2 , and 





1 'an 


Ql2 


Ol3^ 




c ll • 


r (i) \ 


QiQ\ = 


0*12 


«22 


«23 


, Qi = 


c (1) 

c 21 • 


• C 2a« 




\ a 13 


fl 23 


a 33 J 




c (1) 


• C 3 Q W/ 



(27) 



Q2Q2 



^44 045 046^ 

a 45 a 55 °56 
y a 46 a 56 a 66 y 



Q 2 



/ C (2 A 

c (2) 

c 21 

(2) 



(28) 



Here a t = ( N 2 zf), 



0,1,2. Partial transpose on the third qubit of pa results in 
n n n -,,„ n n \ 



a 00 











«12 


«13 











a u 


























«22 


«23 











Cf46 








fl 23 


033 
























a 44 


a 45 



















a* 5 


a 55 


























^66 











a 46 


a 46 















(29) 



In this case, while the density matrix has one zero eigenvalue, the partial transpose has no 
zero eigenvalue. Apart from the eigenvalues an and a 66 , the other six eigenvalues are those 
of the matrices A and B, where: 



^ a 00 a l2 a l3^ 

0^2 ^44 ^45 



a l3 a 45 



, B = 



^22 «23 a 46^ 
a 23 a 33 «56 



055 y 



a 46 a 56 







(30) 
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The characteristic equation of matrix A is 

A 3 — A 2 (a o + ^44 + O55) — -M a 45 + a 12 + a 13 — a 00 a 44 — ^00 a 55 — 0-44^55) (31) 
+ (^00^45 + + 0-44^13 — G[55«00«44 — 2012013045) = 0. 

Here we see that average of coefficients of A 2 and A goes as —1 — 2/N 2 and 4/N 2 while that 
of constant term goes as 16/iV 5 . Thus typically the determinants of A and B, which are the 
negative of the constant term, are negative, and hence the three eigenvalues of each matrix 
can either all be negative or have one lone negative value. That the latter is the case follows 
on noting that the traces of these two matrices are positive. 

One may estimate the negative eigenvalue now. Assuming that the terms containing A 3 
and A 2 are of smaller order, and keeping only the linear and constant terms one gets that 
the negative eigenvalues of A is approximately —4/N 3 . It is immediately verified that the 
assumptions just made are justified. The characteristic equation of matrix B is 

A 3 - A 2 (a 22 + 033) - A(a 2 6 + a 2 6 + a 2 3 - 022^33) + ^22^56 + a 33fl46 _ 2a 23a46a 5 6 = 0. 

The average of the coefficients of A 2 and A go as —2/N and 4/N 2 while that of the constant 
term goes as 16/iV 5 . Using an argument similar to that used to approximate the negative 
eigenvalue of the matrix A, we can approximate the same for matrix B as —4/N 3 . The 
log-negativity is given by E LN w ~ 8/^ 3 + 8/iV 3 = 16/iV 3 . This estimate is 

again in very good agreement with numerical results as shown in Fig. @, including both 
the exponent and the constant. 

V. DENSITY OF STATES BEFORE AND AFTER PT 

The spectral properties of the reduced density matrix which represents the state of the 
block whose entanglement is under investigation is of natural interest. Apart from being 
positive semi-definite the eigenvalues of the reduced density matrix of a pure random state 
has a characteristic distribution or "density of states" , which is discussed further below. In 
contrast the corresponding spectrum for the partial transpose need not be positive semi- 
definite; indeed if the density of states now of the PT of the reduced density matrix, has 
support in the negative numbers, the corresponding state is entangled or NPT. The mecha- 
nism that is responsible for the transitions pointed to above remains to be fully investigated, 
however the density of states of the reduced density matrix may be playing a crucial role. 
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FIG. 5. Density of states of the reduced density matrix pa (left column) and its partial transpose 
p T A (right column). The block length (to) is equal to 6, the total number of qubits (N) is 22 and 
the particle number (I) varies as shown. The insets shows an enlarged view of the region near the 
origin of the respective figures. In the insets of the right column a vertical line at the origin has 
been shown to draw attention to the negative part of the spectrum. 
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To elucidate this, as discussed in the Introduction, when I < m, that is the number 
of particles is smaller than the block length, there are many exact zero eigenvalues in the 
density matrix p A , in contrast when I > m the density matrix becomes of full rank. In fact 
when I = m a detailed study of the density matrix shows a density of states that still diverges 
at 0, while for I > m the density of states vanishes at zero. Negative eigenvalues develop in 
the partial transpose of the rank-deficient matrices corresponding to the case I < m. When 
the density of states of pA is bounded away from zero, as in the case of I > m, the partial 
transpose is also bounded away from zero and has typically only positive eigenvalues. In the 
marginal case when I = m the divergent density of states of pa seems to lead to negative 
partial transpose. Thus whenever a density matrix has exact zero eigenvalues, or has a 
divergent density of states at zero, it will be typically NPT, and hence entangled. As the 
number of particles is increased beyond the block size, the spectrum of pa gets bounded away 
from zero and it becomes PPT. Thus the transition seems to originate in the transition of 
the density of states of the reduced density matrix, which in turn is due to the rank of the 
density matrix becoming full at the point of transition. However we emphasize that the 
observations made here are partially numerical and further work on the partial transpose of 
rank-deficient matrices is necessary to justify them rigorously. 

If \ip) is a full random state of N qubits, mixing all particle numbers together, and let 
the subset A have Hilbert space dimension d,A and the complementary set, dimension d% 
(6,a < ds). Then the density of states of the reduced density matrix of a subset A of the 
cubits pa, if d>A, ds 3> 1, will typically be distributed according to the Marcenko-Pastur rule 
29[: 

&aQ y/ (^max — A) (A — X m in) 



T ( i 2 A ^ < 32 > 



fW 

x = = t a { 1 + q ± 7qP « = *>^- 

In the "symmetric" case of Q = 1, the density of states diverges at the origin, else it is 
bounded away from zero. In fact much is known about the distribution of the smallest 
eigenvalue in the symmetric case, including its distribution and average (l/d A ) [30]. 

Corresponding questions for definite particle subspaces are of natural interest, and we 
present some results here for the density of states, but only in so far as they pertain to the 
problem of entanglement transition studied above. Thus in addition to the density of states, 
-P(A) of the density matrix pA the density of states, Pr(fi), of the partial transpose, p T A is of 
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interest. In Fig. fl5]) is shown the density of states before and after the partial transposition 
for a case when there are N = 22 qubits in all. The block A consists of m = 6 qubits and the 
density of states is shown as the particle number is changed across this value. Looking at 
the density of states -P(A) for the case / = 6 particles one can see the divergence at the origin 
as well as several clumps of eigenvalues. The origin of the clumps is quite easily understood 
as arising from the individual Gk blocks acting as practically independent density matrices, 
the trace normalization condition being the only constraint amongst them. These individual 
blocks then tend to have density of states that are of the nature of the Marcenko-Pastur 
distribution with suitable dimensions. Thus roughly, especially for large N, the density of 
states is pretty much a superposition of such distributions. 

The eigenvalues A of the extreme nonzero blocks Go and G r where r = m or I depending 
on whether I > m or < m, are special in the sense that there is only a lone nonzero eigenvalue 
and they do not follow the Marcenko-Pastur distribution. In any case the "block" G is just 
a number, while G r is a matrix for / < m and a number for / > m. For example in the case 
when m = 2, these two numbers are the values of aoo and 033 of Eq. (|6]), which are seen to 
be the sum of squares of the normally distributed coefficients. Therefore it is easy to see 
that in general they are chi-squared distributed with number of degrees of freedom d: 

2d/2T ! (rf/2) z^-V*/ 2 , x e [0, 00) and d > 1, (33) 

where x = AJ\f. The number of degrees of freedom d depends on whether I < m or / > m. In 
either case for the eigenvalue of Go, d is equal to ( N ~i m ) ■ In the case of the eigenvalue of the 
block G r ,d= (fj™) when I > m, and d = (7) for / < m. Thus when I = m, the number of 
particles is the block size, d = 1 for the eigenvalue of the block G r . Note that when d = 1 
the chi-squared distribution diverges at the origin, unlike the case d > 1. Thus although at 
the transition point / = m, the density matrix is of full rank, it has a divergent density 
of states arising from this eigenvalue. Note that when / > m these lone eigenvalues can 
never lead to a divergent density of states. Also from the inequality in Eq. (J5]) it follows 
that there will never be a symmetric case of the Marcenko-Pastur distribution when / > m. 
Thus indeed this completes the proof that the density of states does not diverge at zero 
when / > m, while it does for / < in. 

Not much is known of the density of states of the partial transpose even for the case o 



full random states, except for a recent mathematical study [3JJ and an ongoing work 32 1 
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which shows for instance that when a density matrix has a symmetric Marcenko-Pastur 
distribution (Q = 1), its partial transpose has a semi-circle distribution. Indeed the density 
of states of the partial transposed matrix p T A as shown in Fig. (J5j) is not very different from 
that of the density matrix itself. The important exceptions are cases where the density of 
states diverges at the origin (in the case of / = 5 and 6) and the density of states of p r A 
clearly has support in the negative numbers, indicating the NPT nature of pa- In contrast 
when there are 7 particles the PT has almost no negative eigenvalues. A much more detailed 
study of the tails of these distributions show a very small fraction of negative eigenvalues, 
indicating that entanglement when present is very rare. This is reflected in the exponentially 
small probability of entangled states after the transition (I > m). Thus the origin of the 
entanglement transitions seems to lie in the change of character of the density of states of 
the reduced density matrix around zero. 



VI. DISCUSSIONS AND CONCLUSION 

While we have focused on the study entanglement within a block of qubits, one may also 
ask if transitions are seen in the entanglement of the block with the rest of the qubits, say 
measured by the von Neumann entropy. Preliminary work not presented here, as well as from 
the discussions above we believe that no such transition is seen. This is due to the fact that 
von Neumann entropy is fairly insensitive to the nature of the density of states around the 
zero eigenvalue. For instance even in the case of full random states, the symmetric states 
(equal bipartitions) do not possess qualitatively different entanglement entropy from the 



non-symmetric ones [33j. On the other hand, entanglement within the block, as measured 
by the log-negativity (or the concurrence in the case of two qubits) is sensitive to the presence 
of a large number of zero or near zero eigenvalues. 

In summary this paper has given definitive evidence of a transition in entanglement 
between two qubits as the number of particles is increased to three. Using log-negativity 
it is shown that the following generalization would hold: the entanglement content in m 
qubits decays algebraically with N, the number of qubits, if the number of particles I < m, 
and exponentially if I > m. Various exponents in the case of algebraic decay have been 
analytically derived for the case of concurrence as well as the log-negativity. The observation 
of a transition is further strengthened by studying the density of states of the reduced density 
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matrix and its partial transpose. It is shown that the rank of the density matrix is not full 
till the number of particles is precisely equal to the block size; and that even at exactly the 
marginal case, the density of states of the reduced density matrix diverges at zero, although 
it is of full rank. The exact zero eigenvalues or the large number of very small ones seem 
to translate on partial transpose to negative eigenvalues, thus leading to typically entangled 
states. This is the case as long as the number of particles is less than or equal to the block 
size. Otherwise the density of states vanishes at zero and leads to a predominantly positive 
partial transpose, which results in the exponentially small entanglement. 

The question of whether the transition studied is also observed on using other entangle- 
ment measures is a natural and interesting one. It is quite easy to see that the negativity 
measure (rather than the log-negativity studied here) also undergoes such a transition. Note 
that unlike the log- negativity measure which is not convex but is nevertheless an entangle- 
ment monotone [34J, the negativity measure is both convex and an entanglement monotone 
22]. While many other measures, such as distillable entanglement [lj], are difficult to com- 



pute, it seems very plausible that the transition is indeed independent of the particular 
measure used. This is strengthened by the study of the spectra of the reduced density ma- 
trix and its partial transpose, and the fact that the transition may well have its origins in 
the behavior of their density of states. 
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